On the convergence of entropy consistent schemes for lubrication type equations in multiple space dimensions

We present nonnegativity-preserving finite element schemes for a general class of thin film equations in multiple space dimensions. The equations are fourth order degenerate parabolic, and may contain singular terms of second order which are to model van der Waals interactions. A subtle discretization of the arising nonlinearities allows us to prove discrete counterparts of the essential estimates found in the continuous setting. By use of the entropy estimate, strong convergence results for discrete solutions are obtained. In particular, the limit of discrete fluxes Mh(Uh)∇Ph will be identified with the fluxM(u)∇(W ′(u)−∆u) in the continuous setting. As a by-product, first results on existence and positivity almost everywhere of solutions to equations with singular lower order terms can be established in the continuous setting.